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Abstract 

(D _ 

I We sliow liow to write a set of brackets for tlie Langevin equation, describing the 

dissipative motion of a classical particle, subject to external random forces. The method 
does not rely on an action principle, and is based solely on the phcnomenological description 
^ , of the dissipative dynamics as given by the Langevin equation. The general expression for 

I the brackets satisfied by the coordinates, as well as by the external random forces, at 

different times, is determined, and it turns out that they all satisfy the Jacobi identity. 



0^ ' Upon quantization, these classical brackets are found to coincide with the commutation 

o 

, to microscopic conservative quantum models for the friction mechanism. 



1 Introduction 

The study of dissipative systems is of fundamental interest in many fields of physics, ranging 
^ ■ from statistical mechanics to condensed matter, atomic physics etc. 

H : 

It is well known that, classically, the action of a heath bath on a particle can be described 
(say, in one space dimension) by the " Langevin force" : 

t 

dt' ii{t-t')±{t') +F{t) , (1.1) 

-co 

where ^{t — t') is a friction coefficient, called "memory function", and F[t) is a random force. 
The memory function is a phcnomenological quantity, depending on the detailed features of the 
coupling with the bath. Apart from the obvious requirement that it should be non-negative at 
all times, it is possible to show 1, that the second law of thermodynamics requires its Fourier 
transform ii{z) to have a real positive part on the real axis: 

Re [/i(Lj + iO+)] > -oo<u;<oo. (1.2) 



*Talk given by the author. 
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The random force F{t) is usually assumed to have zero mean: 

< F{t) >= , 

and to satisfy the Gaussian property, according to which all odd correlators of vanish, while 
the even ones can all be written as sums of products of the two-point function < F(t)F(t') >. 
The expression of the latter function, is determined by the fluctuation-dissipation theorem, in 
terms of the temperature T of the bath, and of the memory function. For example, in the case 
of a memory function of the form /x(t — t') = f6{t — t') (which corresponds to the original form 
of the Langevin equation), one has: 

<F{t)F{t')>=2kTf6{t-t'), (1.3) 

where k is Boltzmann constant. 

If the particle of mass m is subject also to an external conservative force, with potential 
V{x), its motion is then described by the (generalized) "Langevin equation": 

mx + f dt'n{t - t') x{t') + V'{x) = F{t) . (1.4) 

J —oo 

It is well known that this equation describes the approach to equilibrium of the particle 
(if one assumes of course that the heath bath is infinite in size and remains in equilibrium 
at all times), and indeed one can prove that, for any choice of the initial conditions, the 
probability density in the particle's phase space approaches, for large times, the canonical 
Maxwell-Boltzmann distribution. 

It is natural to ask if an analogous theory of dissipation exists, for a quantum particle 
in interaction with a quantum-mechanical bath. In particular, is there a quantum Langevin 
equation, that one can use to describe the influence of the bath on the quantum behavior of 
the particle? At a first glance, it is not very clear what one should expect. Since the Langevin 
force describes the influence of the bath on the particle, and since the bath is itself a quantum 
system, it is quite possible that both the random force F{t) and the memory function /i(t — t') 
become operators, in the quantum theory. The simplest possibility, which we shall assume, is 
to keep /i(t — t') a c-numbcr. However, this is unlikely to be the case for F{f), since the l.h.s. 
of the Langevin equation, which involves the particle's coordinate and momentum, is a-priori 
an operator. In the Quantum theory, besides the randomness already present in the classical 
theory (for T > 0), there will exist a further source of randomness due to intrinsic quantum 
fluctuations. The latter should be encoded by a set of commutators: 

[F{t),F{t')] , [xitinn , m),F{t')]- (1.5) 

Consideration of unequal-time commutators appears to be necessary, because the random force 
F{t) does not obey any equation of motion, and so knowledge of equal-time commutators is 
not sufficient. 
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2 A microscopic model for dissipation 



The standard approach to dissipation in Quantum Mechanics, is based on the physical picture 
that dissipation arises from coupling of the system of interest with a thermal bath. According 
to this picture, one considers a conservative microscopic model for the bath, usually consisting 
of a large number of oscillators, and postulates a certain form for the system-bath interaction. 
Elimination of the degrees of freedom of the bath gives rise to an effective equation of motion for 
the particle, including both the damping force, and the fluctuating force. An implementation 
of this philosophy, in the path-integral formalism, was developed long ago by Feynman and 
Vernon Pj. 

In a canonical framework, Ford et al. rP proposed the following simple independent oscil- 
lators model for the particle-bath system: 



2m ^ > ^ 



L 

One has the standard commutation rules: 



(2.6) 



= ifi , \qj,Vh\ = • (2.7) 

The equations of motion implied by the above hamiltonian are: 

mx + V'(x) = mj ujjiqj - x) (2.8) 
j 

iij + uj]qj = uj]x . (2.9) 

The general solution of Ea. ()2.9() is: 

ft 

qj{t) = q^{t) + x{t) - / dt' cos[ujj{t - t')]±{t') , (2.10) 



where q^(t) is the general solution of the homogeneous equation: 



qj{t) = qj cos{uJjt) + pj . (2.11) 

mjujj 

The choice of the retarded solution in Ea. (|2.1Uj) explicitly breaks time-reversal. The picture 
one has in mind is that the particle is held fixed at x = at t = — oo. Upon plugging Eq. (|2.10|) 
into Ea. H2.8|) . we see that the effective equation of motion for the particle of interest becomes 
a Langevin Equation, with a memory function and a fluctuating force, expressed in terms of 
the oscillator parameters. In detail, one finds: 

/j{t) = mj uj'j cos{LOjt)e{t) , (2.12) 
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where &{t) is the Heaviside function. It is shown in ^ that, by a suitable choice of parame- 
ters, one can reproduce the most general form of the memory function, compatible with the 
principles of thermodynamics. As for the random force, the model gives: 



F{t)=J2m,Lo]q^{t), (2.13) 



and accordingly one finds: 



^[Fit),Fit')] = -J2^, sin[^, (t - t')] = ^{t - t') + ^{f - t) . (2.14) 

j 

The question arises: how much of these commutation relations depend on the detailed micro- 
scopic model? A related question is whether it is possible to obtain the commutators for x{t), 
p{t) and F{t) directly from the Langevin Equation, without making recourse to microscopic 
models? 



3 PB for dissipative systems 

Our approach to the quantization of the Langevin's equation proceeds in a way similar to that 
followed for conservative systems: 

• Step 1: define first a set of Poisson Brackets (PB) { , } on the classical "phase space" V; 

• Step 2: replace PB's by commutators: 

{A,B}^^[A,B] . 
in 

The problems that arise in step 1, are twofold: 

• Problem 1: what is the particle's phase space "P? Since the particle is acted on by a 
random force, initial data do not determine its future evolution. In fact, since F{t) can 
a priori be anything, any path x{t) whatsoever is a possible trajectory for the particle. 
Thus, if we think of the phase space as the set of actual motions of the particle, we need 
to take for V the infinite dimensional set of all paths x{t): 

• Problem 2: dissipative equations cannot be derived, in general, from an action principle. 
How does one get then a PB on 

In order to address both questions, we have elaborated JH' a definition of PB's, that relies 
directly and exclusively on Langevin's equation. To see how this can be done, let us look at the 
conservative limit, when fi(t) = F{t) = 0: 

mi + V'{x) = . (3.15) 
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Then, a canonical Poisson bracket exists, and we can use it to define a retarded two-point 
function: 

G-{t,t') := {x{t),x{t')} ,for t>t' , 

G-{t,t') = , for t < t' , (3.16) 

We can easily derive the equation satisfied by G~{t,t'), for t > t' . Indeed, linearity of the PB, 
and the Leibnitz rule imply: 

"'dt2^"(*'^') = "i^{^(0,3;(t')} = m{x{t),x{t')} = 

= -{V'{x{t)),x{t')} = -V"{xmx{t),x{t')} = -V"{x{t))G-{t,t') . (3.17) 
In view of the boundary condition: 

\im m^G~{t,t') = {mx{t),x{t)} = -1 , (3.18) 
we can write an equation valid for all times: 

(m^ + V"ix{t))^ G-{t,t') = -6{t,t') , (3.19) 

which shows that G~{t,t') is the retarded Green's function for small perturbations of the 
classical motion x{t). One can similarly define an advanced two-point function G~^{t,t'): 

G+{t,t') := -{x{t),x{t')} ,for t < t' , 

G+{t,t') = , for t>t' . (3.20) 

It is easy to verify that G~^{t,t') is the advanced Green's functions for small perturbations of 
classical motions. If we define the "commutator function" G{t,t'): 

G{t, t') = G-{t, t') - G+{t, t') , (3.21) 

we obviously have: 

{x{t),x{t')} = G{t,t') y t,t'. (3.22) 
We see that antisymmetry of the PB follows from the reciprocity relation: 

G+{t,t')=G'{t',t) . (3.23) 

Our proposal |3] to define the PB for the Langevin equation, is to take Ea. ()3.22|) as the 
definition of the PB. This way of defining PB, in terms of the Green's functions associated 
with the operator that describes small perturbations of classical motions, was introduced, for 
systems admitting an action principle, by Peierls |^ (For a review of Peierls brackets, see Refs. 
[HI El-) • Here, we are trying to extend this method to non-conservative systems, for which an 
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action principle cannot be found in general. So, we consider the retarded Green's function for 
the Langevin equation: 

m^ + V"{x{t))^G-{t,t') + j' ^dT^i{t-T)^G-{T,t') = -5{t-t') . (3.24) 

Notice that the lower extremum of integration in the driving term is i', because of the retarded 
boundary condition. Notice also that the "background" curve x{t) is an arbitrary path, since 
any path is a possible motion for the particle, the random force being arbitrary. 

We take G'^(t,t') as the Green's function of the adjoint of Ea. H3.24() : 

+ V"{x{t))] G+{t,t') - f' dTfi{T-t)^G+{T,t') = -6{t-t') . (3.25) 
dt'' J Jt dr 

With this choice, the reciprocity relation Ea. H3.23() still holds. This is easily seen if we write 
Eos. ()3. 241) and (|3.25|) in the form of integral equations (with singular kernels): 

/oo 
dT Lit, t)G- (t, t') = -5{t - t') (3.26) 
-oo 

/oo 
dT L^{t, t)G+{t, t') = -5{t - t') , (3.27) 
-oo 

where the superscript T denotes transpose (transpose is the same as adjoint, because we are 
in the real field). We switch to DeWitt |5j condensed notation, in which the continuous time 
variable is treated as a discrete index and so linear integro-differential operators are written 
as matrices. In this notation: 

LijG-^^ = -5^i , (L^)i,G+^'= = . (3.28) 

Multiplication of the second equation by gives: 

G-'\L^)ijG^^^ = G''^6^ = G~^^ . (3.29) 

However, using the first of Ea. (|3.28|) . the l.h.s. is also equal to: 

G-'\L^)ijG+^^ = G-'^LjiG+i^ = 5\G+^^ = . (3.30) 

Comparison of the r.h.s. of the above two Eqs. proves the reciprocity relation. 
The important issue is to check the Jacobi identity. A direct computation gives: 

x^}, x''} + cycl.perm = G'^&^j + &G''ii + G^^G'\i = T'^^ , (3.31) 

and so the Jacobi identity is fulfilled is the quantity T*-'^ vanishes. Use of reciprocity relations 
allows to write T*-'^ only in terms of G^*-' , and its functional derivatives. The latter are easily 
obtained from Ea. (|05)) : 

Uj^iG-^^ + U,G-^\i = Q. (3.32) 
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Multiplication by G^™, and use of reciprocity relation, gives: 

^-mfe^ ^ _^-im^.^.^^-jfc ^ (3.33) 

A simple computation then gives: 

T^^ = (G-'*G-"^^G-"^ + c.p.){L^n - Km), I (3.34) 

This shows that the Jacobi identity holds if and only if the antisymmetric part of the operator 
for the perturbations, L, has vanishing functional derivatives (namely is independent on the 
unperturbed path x{t)). This is surely the case if the equations of motion derive from an action, 
because then Lnm is the second functional derivative of the action S, and since derivatives 
commute, it is symmetric: 

Lmn — S^rnn — S^nm — Lmn ■ (3.35) 

Now, the operator L for the Langevin equation is not symmetric, however its antisymmetric 
part has vanishing functional derivatives, because of the linear character of the friction term. 
Indeed, if we restore the plain notation for integral operators, we can write the antisymmetric 
part of L (acting on a function ip{t)) as: 

{L - L^)V'(i) = / dt' [fi{t - t') + fi{t' - t)]i!{t') (3.36) 



The kernel of L — is independent on the unperturbed path x{i) and so it has null functional 
derivatives. Then, the friction term does not spoil the Jacobi identity. 

By taking the Langevin Equation as a definition of F{t), we can evaluate the PB satisfied by 
the external force: 

{F(t),F(0} = ^(t-0 + f (t'-t). (3.37) 

This is the same as {l/{ifi)) times the commutator that was obtained in by using an 
explicit microscopic model. It is also possible to verify that the equal-time PB of the particle's 
coordinate and velocity remain canonical: 

{x(t),x(t)} = {i(t),x(t)} = 0, 

m{x{t),x{t)} = 1 . 

It is important to observe that, when friction is present, F{t) has non trivial brackets, and 
so it is inconsistent to set F{t) = 0. This means that a consistent PB can be written for the 
Langevin equation, only if the particle is acted on by an external force. However, if the friction 
term is zero, F{t) has vanishing brackets with everything (including the particle's coordinate); 
then, it is possible to take F{t) = and to restrict the bracket onto the set of solutions of the 
equations of motion. In this way we recover back the (finite dimensional) phase space of the 
conservative system, and its canonical PB. 
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4 Conclusions 



We have defined a set of PB for the generahzed Langevin equation. Distinctive features of our 
approach are: 

• the method does not require an action principle, and is based directly on the equations 
of motion; 

• it can be applied to dissipative equations, for which an action principle does not exist, like 
the Langevin equation. In this way, we obtained the PB directly from the macroscopic 
description of dissipation, as provided by the Langevin equation, without making recourse 
to microscopic models; 

• when dissipation is present, the relevant phase space is, a priori, the space of all paths, 
which is an infinite-dimensional space; 

• in the absence of dissipation, when the system is conservative, it is possible to set the 
random force to zero and restrict the bracket onto the standard phase space, spanned by 
classical solutions of the e.o.m.. One recovers then the usual canonical PB. 

G.B., G.E. and G.M. acknowledge partial financial support by PRIN 2002 SINTESI. 
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